Genetic diversity among the present Japanese population: evidence from genotyping of human cell lines established in Japan.

Japan is often assumed to have a highly homogeneous ethnic population, because it is an island country. This is evident in human cell lines collected from cell banks; however, these genotypes have not been thoroughly characterized. To examine the population genotypes of human cell lines established in Japan, we conducted SNP genotyping on 57 noncancerous cell lines and 43 lung cancer cell lines. Analysis of biogeographic ancestry revealed that 58 cell lines had non-admixed Japanese genotypes, 21 cell lines had an admixture of Japanese and East Asian genotypes, and the remaining 21 cell lines had East Asian genotypes. The proportion of non-admixed Japanese genotypes was similar between lung cancer and noncancerous cell lines, suggesting that patients in Japan may not exclusively have Japanese genotypes. This could influence the incidence of inherited diseases and should be taken into account in personalized medicine tailored to genetic background. The genetic makeup of the present-day Japanese population cannot be fully explained by the ancestral Jomon and Yayoi lineages. Instead, it is necessary to consider a certain level of genetic admixture between Japanese and neighboring Asian populations. Our study revealed genetic variation among human cell lines derived from Japanese individuals, reflecting the diversity present within the Japanese population.

The Cytogenetic Map of the Nile Crocodile (Crocodylus niloticus, Crocodylidae, Reptilia) with Fluorescence In Situ Localization of Major Repetitive DNAs.

Tandemly arranged and dispersed repetitive DNA sequences are important structural and functional elements that make up a significant portion of vertebrate genomes. Using high throughput, low coverage whole genome sequencing followed by bioinformatics analysis, we have identified seven major tandem repetitive DNAs and two fragments of LTR retrotransposons in the genome of the Nile crocodile (Crocodylus niloticus, 2n = 32). The repeats showed great variability in structure, genomic organization, and chromosomal distribution as revealed by fluorescence in situ hybridization (FISH). We found that centromeric and pericentromeric heterochromatin of C. niloticus is composed of previously described in Crocodylus siamensis CSI-HindIII and CSI-DraI repetitive sequence families, a satellite revealed in Crocodylus porosus, and additionally contains at least three previously unannotated tandem repeats. Both LTR sequences identified here belong to the ERV1 family of endogenous retroviruses. Each pericentromeric region was characterized by a diverse set of repeats, with the exception of chromosome pair 4, in which we found only one type of satellite. Only a few repeats showed non-centromeric signals in addition to their centromeric localization. Mapping of 18S-28S ribosomal RNA genes and telomeric sequences (TTAGGG)n did not demonstrate any co-localization of these sequences with revealed centromeric and pericentromeric heterochromatic blocks.

Observations on chromosome-specific sequencing for the construction of cross-species chromosome homology maps and its resolution of human:alpaca homology.

BACKGROUND: The history of comparative chromosome mapping is briefly reviewed, with discussion about the problem that occurs in chromosome painting when size heteromorphisms between homologues cause contamination in chromosomes sorted by flow cytometry that are used in the preparation of chromosome-specific DNA probes. MAIN BODY: As an example, we show in the alpaca (Vicagna pacos) that sequencing of contaminated chromosome sorts can reveal chromosome homologies from alignment with human and mouse genome reference sequences. The procedure identifies syntenic blocks of DNA separated in the human karyotype that are associated in the closely related alpaca and dromedary (Camelus dromedarius) karyotypes. This example provides proof of principal for the validity of the method for comparative chromosome mapping. CONCLUSION: It is suggested that the approach presented here may have application in the construction of comparative chromosome maps between distantly related taxa, such as monotremes and mammals.

THP-1 reference data: Proposal of an in vitro branched evolution model for cancer cell lines.

THP-1 is a representative leukemia cell line and is registered with four different numbers in JCRB and RIKEN BRC cell banks. However, differences between these four lines remain unclear. In our study, these four THP-1 cell lines, JCRB0112, JCRB0112.1 (corresponding to ATCC TIB-202), RCB1189 (DSMZ ACC-16) and RCB3686, have been compared at chromosome and DNA sequence levels. Our results reveal that ploidy has been changed in JCRB0112 and RCB1189, which are triploid and tetraploid, respectively. Patterns of variant frequencies from target sequencing are unique to each ploidy, estimating whole genomic status based on partial sequence data. SNP microarrays showed four distinct profiles with a large-scale loss of heterozygosity, reflected in subtle differences in STR genotypes. Transcriptome patterns suggest that JCRB0112.1 has diverged highly from the other three lines. RCB1189 and JCRB0112.1 responded to PMA faster than RCB3686 and JCRB0112. We have identified RCB3686 as the closest to the original THP-1, which can be an optimal model of AML-M5. These four THP-1 genomes and transcriptomes exhibit significant differences, indicating four independent sublines and demonstrating the influence of genetic drift on gene expression. As these cells share the same name, THP-1 must be accompanied by their registration number of each cell repository. Our data provide genomic features of four THP-1 sublines and serve as a reference profile to classify widely spread THP-1 progenies, which could be distinguished by a comparison of 24 STR markers. Multiple sublines can be generated by separate cell cultures, which would be explained by in vitro branched evolution.

Whole-Genome Sequencing of Vero E6 (VERO C1008) and Comparative Analysis of Four Vero Cell Sublines.

The Vero cell line is an immortalized cell line established from kidney epithelial cells of the African green monkey. A variety of Vero sublines have been developed and can be classified into four major cell lineages. In this study, we determined the whole-genome sequence of Vero E6 (VERO C1008), which is one of the most widely used cell lines for the proliferation and isolation of severe acute respiratory syndrome coronaviruses (SARS-CoVs), and performed comparative analysis among Vero JCRB0111, Vero CCL-81, Vero 76, and Vero E6. Analysis of the copy number changes and loss of heterozygosity revealed that these four sublines share a large deletion and loss of heterozygosity on chromosome 12, which harbors type I interferon and CDKN2 gene clusters. We identified a substantial number of genetic differences among the sublines including single nucleotide variants, indels, and copy number variations. The spectrum of single nucleotide variants indicated a close genetic relationship between Vero JCRB0111 and Vero CCL-81, and between Vero 76 and Vero E6, and a considerable genetic gap between the former two and the latter two lines. In contrast, we confirmed the pattern of genomic integration sites of simian endogenous retroviral sequences, which was consistent among the sublines. We identified subline-specific/enriched loss of function and missense variants, which potentially contribute to the differences in response to viral infection among the Vero sublines. In particular, we identified four genes (IL1RAP, TRIM25, RB1CC1, and ATG2A) that contained missense variants specific or enriched in Vero E6. In addition, we found that V739I variants of ACE2, which functions as the receptor for SARS-CoVs, were heterozygous in Vero JCRB0111, Vero CCL-81, and Vero 76; however, Vero E6 harbored only the allele with isoleucine, resulting from the loss of one of the X chromosomes.

TK6 genome profile compared with WIL2-NS: Reference data to improve the reproducibility of genotoxicity studies.

The TK6 cell line has been widely used for genotoxicity screening by taking advantage of the TK1 mutations. A number of publications have employed TK6 as a tool for the assay; however, its characterization is limited. Because genomes of cell lines are often changed during cell culture, differences in TK6 genomes could be observed between laboratories, which would cause potential problems with reproducibility using TK6. In this study, TK6(IVGT) (JCRB1435), registered with the JCRB Cell Bank as the standard for TK6, has been characterized by karyotyping, SNP microarray and targeted sequencing, and were then compared to WIL2-NS (JCRB9063), a subline derived from a common ancestor with TK6. Changes at the chromosome level were quantitatively assessed by the microarray data, which were shown by DNA sizes and Scales of Genomic Alterations (SGA). An extinct common ancestral cell line, WI-L2 originated from hereditary spherocytosis, has been extrapolated from analysis of genomic signatures shared between TK6 and WIL2-NS, revealing a point mutation in SPTA1, related to the disease. DNA size is calculated to be 102.6 %, 103.1 % and 103.9 % for WI-L2, TK6(IVGT) and WIL2-NS respectively, with the reference values of a normal diploid genome, and each genome shows SGA of 2.8 %, 4.5 % and 4.2 % respectively. Differences between TK6(IVGT) and WIL2-NS are due to regions and sizes of gains, implying that genomes of both cells tend to increase their DNA size independently. Our data provide fundamental genomic features of TK6 and serve as a reference profile of the standardized TK6 cell line, which leads to an increase in robustness of assays using the TK6 cells.

Kasumi leukemia cell lines: characterization of tumor genomes with ethnic origin and scales of genomic alterations.

Kasumi-1 has played an important role in an experimental model with t(8;21) translocation, which is a representative example of leukemia cell lines. However, previous studies using Kasumi-1 show discrepancies in the genome profile. The wide use of leukemia cell lines is limited to lines that are well-characterized. The use of additional cell lines extends research to various types of leukemia, and to further explore leukemia pathogenesis, which can be achieved by uncovering the fundamental features of each cell line with accurate data. In this study, ten Kasumi cell lines established in Japan, including five that were previously unknown, have been characterized by SNP microarray and targeted sequencing. SNP genotyping suggested that the genetic ancestry in four of the ten Kasumi cell lines was not classified as Japanese but covered several different east-Asian ethnicities, suggesting that patients in Japan are genetically diverse. TP53 mutations were detected in two cell lines with complex array profiles, indicating chromosomal instability (CIN). A quantitative assessment of tumor genomes at the chromosomal level was newly introduced to reveal total DNA sizes and Scales of Genomic Alterations (SGA) for each cell line. Kasumi-1 and 6 derived from relapsed phases demonstrated high levels of SGA, implying that the level of SGA would reflect on the tumor progression and could serve as an index of CIN. Our results extend the leukemia cellular resources with an additional five cell lines and provide reference genome data with ethnic identities for the ten Kasumi cell lines.

Squamate Chromosome Size and GC Content Assessed by Flow Karyotyping.

Chromosome homologies in reptiles have been investigated extensively by gene mapping and chromosome painting. Relative chromosome size can be estimated roughly from conventional karyotypes, but chromosome GC content cannot be evaluated by any of these approaches. However, GC content can be obtained by whole-genome sequencing, although complete data are available only for a limited number of reptilian species. Chromosomes can be characterized by size and GC content in bivariate flow karyotypes, in which the distribution of peaks represents the differences. We have analysed flow karyotypes from 9 representative squamate species and show chromosome profiles for each species based on the relationship between size and GC content. Our results reveal that the GC content of macrochromosomes is invariable in the 9 species. A higher GC content was found in microchromosomes, similar to profiles previously determined in crocodile, turtle, and chicken. The findings suggest that karyotype evolution in reptiles is characterized by unique features of chromosome GC content.

Marsupial chromosome DNA content and genome size assessed from flow karyotypes: invariable low autosomal GC content.

Extensive chromosome homologies revealed by cross-species chromosome painting between marsupials have suggested a high level of genome conservation during evolution. Surprisingly, it has been reported that marsupial genome sizes vary by more than 1.2 Gb between species. We have shown previously that individual chromosome sizes and GC content can be measured in flow karyotypes, and have applied this method to compare four marsupial species. Chromosome sizes and GC content were calculated for the grey short-tailed opossum (2n = 18), tammar wallaby (2n = 16), Tasmanian devil (2n = 14) and fat-tailed dunnart (2n = 14), resulting in genome sizes of 3.41, 3.31, 3.17 and 3.25 Gb, respectively. The findings under the same conditions allow a comparison between the four species, indicating that the genomes of these four species are 1-8% larger than human. We show that marsupial genomes are characterized by a low GC content invariable between autosomes and distinct from the higher GC content of the marsupial × chromosome.

Screening for 15 pathogenic viruses in human cell lines registered at the JCRB Cell Bank: characterization of in vitro human cells by viral infection.

Human cell lines have been used in a variety of research fields as an in vitro model. These cells are all derived from human tissue samples, thus there is a possibility of virus infection. Virus tests are routinely performed in clinical practice, but are limited in cell lines. In this study, we investigated 15 kinds of viruses in 844 human cell lines registered at the Japanese Collection of Research Bioresources (JCRB) Cell Bank. Our real-time PCR analysis revealed that six viruses, EBV, HTLV-1, HBV, B19V, HHV-6 and HHV-7, were detected in 43 cell lines. Of them, 20 cell lines were transformed by intentional infection in vitro with EBV or HTLV-1. Viruses in the other 23 cell lines and one EBV transformed cell line are derived from an in vivo infection, including five de novo identifications of EBV, B19V or HHV-7 carriers. Among them, 17 cell lines were established from patients diagnosed with virus-associated diseases. However, the other seven cell lines originated from in vivo cells unrelated to disease or cellular tropism. Our approach to screen for a set of 15 viruses in each cell line has worked efficiently to identify these rare cases. Virus tests in cell lines contribute not only to safety assessments but also to investigation of in vivo viral infection which can be a characteristic feature of cell lines.

HuH-7 reference genome profile: complex karyotype composed of massive loss of heterozygosity.

Human cell lines represent a valuable resource as in vitro experimental models. A hepatoma cell line, HuH-7 (JCRB0403), has been used extensively in various research fields and a number of studies using this line have been published continuously since it was established in 1982. However, an accurate genome profile, which can be served as a reliable reference, has not been available. In this study, we performed M-FISH, SNP microarray and amplicon sequencing to characterize the cell line. Single cell analysis of metaphases revealed a high level of heterogeneity with a mode of 60 chromosomes. Cytogenetic results demonstrated chromosome abnormalities involving every chromosome in addition to a massive loss of heterozygosity, which accounts for 55.3% of the genome, consistent with the homozygous variants seen in the sequence analysis. We provide empirical data that the HuH-7 cell line is composed of highly heterogeneous cell populations, suggesting that besides cell line authentication, the quality of cell lines needs to be taken into consideration in the future use of tumor cell lines.

Homologue-specific chromosome sequencing characterizes translocation junctions and permits allelic assignment.

Chromosome translocations can be detected by cytogenetic analysis, but it is hard to characterize the breakpoints at the sequence level. Chromosome sorting by flow cytometry produces flow karyotypes that enable the isolation of abnormal chromosomes and the generation of chromosome-specific DNA. In this study, a derivative chromosome t(9; 14) and its homologous normal chromosomes 9 and 14 from the Ishikawa 3-H-12 cell line were sorted to collect homologue-specific samples. Chromosome sequencing identified the breakpoint junction in the der(9) at 9p24.3 and 14q13.1 and uncovered the formation of a fusion gene, WASH1-NPAS3. Amplicon sequencing targeted for neighbouring genes at the fusion breakpoint revealed that the variant frequencies correlate with the allelic copy number. Sequencing of sorted chromosomes permits the assignment of allelic variants and can lead to the characterization of abnormal chromosomes. We show that allele-specific chromosome sequencing of homologues is a robust technique for distinguishing alleles and this provides an efficient approach for the comprehensive analysis of genomic changes.

Novel endogenous simian retroviral integrations in Vero cells: implications for quality control of a human vaccine cell substrate.

African green monkey (AGM)-derived Vero cells have been utilized to produce various human vaccines. The Vero cell genome harbors a variety of simian endogenous type D retrovirus (SERV) sequences. In this study, a transcriptome analysis showed that DNA hypomethylation released the epigenetic repression of SERVs in Vero cells. Moreover, comparative genomic analysis of three Vero cell sublines and an AGM reference revealed that the genomes of the sublines have ~80 SERV integrations. Among them, ~60 integrations are present within all three cell sublines and absent from the reference sequence. At least several of these integrations consist of complete SERV proviruses. These results strongly suggest that SERVs integrated in the genome of Vero cells did not retrotranspose after the establishment of the cell lineage as far as cells were maintained under standard culture and passage conditions, providing a scientific basis for controlling the quality of pharmaceutical cell substrates and their derived biologics.

A collection of XY female cell lines.

Discordance between sexual phenotype and the 46,XY sex chromosome complement may be found in certain disorders of sexual development (DSD). Many of these DSD patients with female external genitalia and secondary sex characteristics have undescended testes and male internal genitalia. Causative mutations involving genes of the sex determining pathway, including the androgen receptor, SRY and the 5-alpha-reductase genes, are well-known, but the origin of other cases remain unresolved. In this report, we introduce our collection of lymphoblastoid lines derived from female patients with a 46,XY karyotype. These cell lines have been deposited and registered with the JCRB Cell Bank. They are available for comparison with other DSD cases and for further characterization of genetic loci involved in the mammalian sex determining pathway.

The human vascular endothelial cell line HUV-EC-C harbors the integrated HHV-6B genome which remains stable in long term culture.

Human herpes virus 6 (HHV-6) is a common human pathogen that is most often detected in hematopoietic cells. Although human cells harboring chromosomally integrated HHV-6 can be generated in vitro, the availability of such cell lines originating from in vivo tissues is limited. In this study, chromosomally integrated HHV-6B has been identified in a human vascular endothelial cell line, HUV-EC-C (IFO50271), derived from normal umbilical cord tissue. Sequence analysis revealed that the viral genome was similar to the HHV-6B HST strain. FISH analysis using a HHV-6 DNA probe showed one signal in each cell, detected at the distal end of the long arm of chromosome 9. This was consistent with a digital PCR assay, validating one copy of the viral DNA. Because exposure of HUV-EC-C to chemicals did not cause viral reactivation, long term cell culture of HUV-EC-C was carried out to assess the stability of viral integration. The growth rate was altered depending on passage numbers, and morphology also changed during culture. SNP microarray profiles showed some differences between low and high passages, implying that the HUV-EC-C genome had changed during culture. However, no detectable change was observed in chromosome 9, where HHV-6B integration and the viral copy number remained unchanged. Our results suggest that integrated HHV-6B is stable in HUV-EC-C despite genome instability.

Incorrect strain information for mouse cell lines: sequential influence of misidentification on sublines.

Misidentification or cross-contamination of cell lines can cause serious issues. Human cell lines have been authenticated by short tandem repeat profiling; however, mouse cell lines have not been adequately assessed. In this study, mouse cell lines registered with the JCRB cell bank were examined by simple sequence length polymorphism (SSLP) analysis to identify their strains. Based on comparisons with 7 major inbred strains, our results revealed their strains in 80 of 90 cell lines. However, 12 of the 80 cell lines (15%) were found to differ from registered information. Of them, 4 cell lines originated from the same mouse, which had been generated through mating between two different inbred strains. The genotype of the mouse sample had not been examined after the backcross, leading to strain misidentification in those cell lines. Although 8 other cell lines had been established as sublines of a BALB/c cell line, their SSLP profiles are similar to a Swiss cell line. This affects differences in genotypes between inbred and outbred strains. Because the use of inbred samples and interbreeding between strains are not involved in human materials, our results suggest that the cause and influence of misidentification in mouse cell lines are different from those in human.

Changes of heterogeneous cell populations in the Ishikawa cell line during long-term culture: Proposal for an in vitro clonal evolution model of tumor cells.

Genomic changes in tumor cell lines can occur during culture, leading to differences between cell lines carrying the same name. In this study, genome profiles between low and high passages were investigated in the Ishikawa 3-H-12 cell line (JCRB1505). Cells contained between 43 and 46 chromosomes and the modal number changed from 46 to 45 during culture. Cytogenetic analysis revealed that a translocation t(9;14), observed in all metaphases, is a robust marker for this cell line. Single-nucleotide polymorphism microarrays showed a heterogeneous copy number in the early passages and distinct profiles at late passages. These results demonstrate that cell culture can lead to elimination of ancestral clones by sequential selection, resulting in extensive replacement with a novel clone. Our observations on Ishikawa cells in vitro are different from the in vivo heterogeneity in which ancestral clones are often retained during tumor evolution and suggest a model for in vitro clonal evolution.

The genome landscape of the african green monkey kidney-derived vero cell line.

Continuous cell lines that originate from mammalian tissues serve as not only invaluable tools for life sciences, but also important animal cell substrates for the production of various types of biological pharmaceuticals. Vero cells are susceptible to various types of microbes and toxins and have widely contributed to not only microbiology, but also the production of vaccines for human use. We here showed the genome landscape of a Vero cell line, in which 25,877 putative protein-coding genes were identified in the 2.97-Gb genome sequence. A homozygous ∼9-Mb deletion on chromosome 12 caused the loss of the type I interferon gene cluster and cyclin-dependent kinase inhibitor genes in Vero cells. In addition, an ∼59-Mb loss of heterozygosity around this deleted region suggested that the homozygosity of the deletion was established by a large-scale conversion. Moreover, a genomic analysis of Vero cells revealed a female Chlorocebus sabaeus origin and proviral variations of the endogenous simian type D retrovirus. These results revealed the genomic basis for the non-tumourigenic permanent Vero cell lineage susceptible to various pathogens and will be useful for generating new sub-lines and developing new tools in the quality control of Vero cells.

Novel tools for characterising inter and intra chromosomal rearrangements in avian microchromosomes.

Avian genome organisation is characterised, in part, by a set of microchromosomes that are unusually small in size and unusually large in number. Although containing about a quarter of the genome, they contain around half the genes and three quarters of the total chromosome number. Nonetheless, they continue to belie analysis by cytogenetic means. Chromosomal rearrangements play a key role in genome evolution, fertility and genetic disease and thus tools for analysis of the microchromosomes are essential to analyse such phenomena in birds. Here, we report the development of chicken microchromosomal paint pools, generation of pairs of specific microchromosome BAC clones in chicken, and computational tools for in silico comparison of the genomes of microchromosomes. We demonstrate the use of these molecular and computational tools across species, suggesting their use to generate a clear picture of microchromosomal rearrangements between avian species. With increasing numbers of avian genome sequences that are emerging, tools such as these will find great utility in assembling genomes de novo and for asking fundamental questions about genome evolution from a chromosomal perspective.

Afrotheria genome; overestimation of genome size and distinct chromosome GC content revealed by flow karyotyping.

Afrotheria genome size is reported to be over 50% larger than that of human, but we show that this is a gross overestimate. Although genome sequencing in Afrotheria is not complete, extensive homology with human has been revealed by chromosome painting. We provide new data on chromosome size and GC content in four Afrotherian species using flow karyotyping. Genome sizes are 4.13 Gb in aardvark, 4.01 Gb in African elephant, 3.69 Gb in golden mole and 3.31 Gb in manatee, whereas published results show a mean of 5.18 Gb for Afrotheria. Genome GC content shows a negative correlation with size, indicating that this is due to differences in the amount of AT-rich sequences. Low genome GC content and small variance in chromosome GC content are characteristic of aardvark and elephant and may be associated with the high degree of conserved synteny, suggesting that these are features of the Afrotherian ancestral genome.